1 




THE EFFECTS OF EXPLOSIVE SHOCK WAVE PROPAGATION 
THROUGH A SOLID STATE MOLECULAR STRUCTURE 



THE EFFECTS OF EXPLOSIVE SHOCK WAVE PROPAGATION 
THROUGH A SOLID STATE MOLECULAR STRUCTURE 



by 

David Edwin Clark, B.S. 

/ ' t e • 



THESIS 

Presented to the Faculty of the Graduate School of 
The University of Texas at Austin 
in Partial Fulfillment 
of the Requirements 
for the Degree of 

MASTER OF ARTS 

THE UNIVERSITY OF TEXAS AT AUSTIN 
August 1988 



TABLE OF CONTENTS 



Chapter Page 

List of Figures iv 

I . Introduction 1 

II. Method 5 

III. Discussion 26 

IV. Conclusions 33 

Tables 36 

Figure Captions 39 

Figures 46 

References 125 

Vita 127 



iii 



LIST OF FIGURES 



Figures 1-2 


Model 


Figures 3-5 


Classical Positions 


Figures 6- 1 1 


Classical Bond Energies 


Figures 12-33 


Quantum State Occupation Probabilities 


Figures 34-61 


Quantum Atom Positions 


Figures 62-69 


Mode Energies 


Figures 70-79 


Quantum Bond Energies 



IV 



I. Introduction 

Shock waves have been studied in the past in order to provide insight to the 
energy released during an external shock to a system and the resultant physical and 
chemical reactions which occur. A detonation may be defined as a supersonic wave 
which propagates throughout a crystal lattice structure, and whose energy is released 
during an exothermic chemical reaction. ^ It is convenient to classify a detonation into 
two categories, the first category being defined as the shock front which is 
characterized by a sharp rise in energy and agitation over time occuring on the order of 
picoseconds. As the shock front passes through a lattice, an almost instantaneous 
change in atom position and momentum may be observed. The second category is the 
shock wave which is usually characterized by the slowly decaying nature of the wave 

9 

until thermal equilibrium has been established.” The shock wave transfers energy 
through the lattice in a much slower fashion. Time delays are experienced with respect 

to atom displacement as well as energy transfer. Shock waves propagate through most 
materials on the order of 2500 meters per second ( 5.49 x 10'“ a^/ps ).^ There has 

been considerable attention in the past devoted to the examination of shock waves and 
whether they may be classified into two seperate categories ( first and second sound ). 
This idea has been disputed in the literamre for several years. 

Shock waves produce a variety of anomolous optical and electrical properties as 
well as non thermal chemistry, and they are capable of producing bond dissociation, 
phase changes, and a redistribution of vibrational, electronic and rotational energy 
within a lattice. Shock waves are also unique in that energy may be deposited to 
vibrational modes of a molecule under nonequilibrium conditions. ^ 
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Early studies concerning simulations of shock waves propagating within crystal 
lattices were restricted to classical dynamics with finite dimensions, consisting of non- 
oscillatory atoms prior to the introduction of a shock wave.^"^ These studies were 
primarily aimed at understanding the shock wave, and its velocity and kinetic energy 
effects upon a crystal lattice. Later studies^'^ devoted considerable attention to the 
characteristics of the shock wave and the various phases of velocity transformations. 
Classical studies have been completed on systems which were initially at a nonzero 
temperature (zero temperature being defined as no oscilatory motion),9-l0 and on 
models which incorporated deformations within the crystal lattice.l^'^'^ These studies 
were concerned with affects of free radicals, mass defects which were randomly 
dispersed, and mass defects of heavy and light impurities on shock wave propagation 
throughout the lattice. 

There have been few quantum mechanical studies concerning shock waves. 
Dancz and Rice derived expressions for the quantum mechanical equations of motion 
for coupled anharmonic oscillators.^^ TTiere have been other studies which were also 
devoted to deriving equations of motion for quantum systems by the introduction of 
raising and lowering operators for interacting Morse potentials. Other studies 
examined the probability of occupation of excited states and the state to state energy 
flow in which energy became trapped in excited states. 

Experimental work using picosecond lasers has also been done in order to 

understand shock wave propagation through lattices. In this expenment, a 

1 2 2 

picosecond laser ( capable of 10 Watts/cm ) imposed a shock upon a sample of 
water which was analyzed using Raman Spectroscopy. This was the first experiment 
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in which energy transfer from a macroscopic shock to intermolecular states has been 
observed. 

Most of the previous studies have been devoted extensively to the classical study 
of shock waves, and the subsequent interface with a lattice model using classical 
mechanics. Much work has been devoted to examining the shock wave and the 
resultant energy transfer associated with it. Recent studies fWyatt and Marston] have 
been devoted to describing a shock wave effect quantum mechanically in a one 
dimensional crystal lattice. This work was primarily devoted towards a high energy 
shock wave and its effect upon a completely harmonic system. A classical linear model 
of twenty diatoms was formulated and perturbed via a ballistic particle of equivalent 
mass. The inner two diatom pairs were chosen for the classical study in order to 
restrict center of mass movement within the lattice and to represent inner lattice 
substituents. These inner atoms were perturbed by direct interactions of their neighbor 
atoms (entering and absorbing driving atoms). Wyatt and Marston used the action of 
these atoms to establish a time-dependent driving potential for two diatoms. The 
entering atom (q<) provided the initial interaction from the shock to the cluster, 
whereas, the absorbing atom (q>) acted as a resevoir for the energy to be released from 
the cluster. The dynamical data from the classical results of q< and q> was used 
directly to formulate an interaction potential for a quantum mechanical model of two 
diatoms; thus, the classical and quantum model experienced similar potentials. The 
quantum model was represented by four normal modes, upon perturbation yielded 
quantum dynamical data for comparison to the classical system. Wyatt and Marston 
showed that a high energy shock imposed upon a harmonic quantum system 
converges, as expected, to the classical results. 



This work is devoted to the replacement of harmonic oscillators with cubic 
anharmonic oscillators and in making comparisons to classical d\Tiamics. A similar 
anharmonic classical model is established in order to provide for a driving potential in 
our quantum system. In the quantum system the cubic anharmonicities provide various 
mode couplings; thus, our Hamiltonian may be not separated exclusively into the four 
modes as in Marston's; however, it wall be represented by all four modes. This fact 
will restrict our quanmm basis size. 

The remainder of this thesis is organized in the following fashion. Section II 
develops the model for both classical and quantum systems. The classical equations of 
motion are derived showing how w'e obtained d>Tiamical data for comparison and for 
the quantum potential. The quantum coupled anharmonic Hamiltonian is derived using 
the normal modes representation and the time dependent system is developed by 
approximately solving the Schrodinger time dependent equation. Quantum formulation 
for position expectation values, normal mode energies, and bond energies is presented 
in this section. 

Section rn presents results of quantum harmonic and anharmonic calculations 
of several basis sizes and compares these results to classical d\'namical data. Section 
rv draws conclusions based upon thse results and also provides commentary' based 
upon different formulations of the problem. 



II. Method 



Model 

A one-dimensional twenty-diatom classical model was established in order to 
study the effects of the propagation of a shock wave through a crystal lattice. All of 
the atoms in the chain are of equal mass (mass of fluorine atom); the intermolecular 
distances represent the intermolecular separations of diatomic fluorine. Each diatom 
was connected with a force constant K> and equilibrium separation po. The 
equilibrium separation between diatoms was R<j with a connecting force constant K< 
(see Figure 1 and Table 1). A ballistic particle of equal mass will impose a shock wave 
upon the chain of atoms. In order to compare a classical model with a quantum model, 
the central pair of diatoms in the classical one dimensional chain will be isolated. This 
pair of diatoms is far enough from the ballistic particle to reduce any sporadic results 
from the strong potential and a smaller system of two diatoms was chosen for 
comparison to an equivalent quantum model based on the fact that a twenty diatom 
cluster is numerically intractable to handle quantum mechanically. Classically, the 
central diatom pair wiU be perturbed via interaction with the nearest atoms in the chain. 
The entering atom (q<) provides the initial interaction from the shock to the cluster; 
whereas, the absorbing atom (q>) acts as a reservoir for the energy to be released from 
the cluster. The dynamical data from the classical results of q<(t) and q>(t) will be used 
directly to formulate an interaction potential for a quantum mechanical model of two 
diatoms; thus, the classical and quantum models will be experiencing similar 
potentials, thus allowing for direct comparisons. 

Classically, these atoms may be thought of as forty harmonic or anharmonic 
oscillators. At time t=0, atom positions were arbitrarily set at their equilibrium values. 
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and atom number one in the chain was originally chosen to have zero momentum. 

Atom momenta are easily obtainable through the following relationships: 

(?2 ■ P i) 1 2 3 

+ iK^(q2-qi -p/ + Y(q2-qi-Po) =Et 
2p 



( 1 ) 



Pi = momentum of atom number one = 0 
po = equilibrium separation (Table 1) 

(q2 ■ qi) = poi since at equilibrium 

Et = Total energy (sum of potential, and kinetic terms) 

|i = reduced mass (Table 1) 

Y = degree of anharmonicity (Table 1) 

K> = Force constant between atoms (Table 1) 

The total energy (Ej) between any two atoms was calculated from the ground state 
energy of a diatom pair (Ej = 2 hco where co =VK/|i.). Equation ( 1 ) readily reduces to 



(P2) 

2p 



- Et 



since the potential energy term and the anharmonic term consistently drop 



out of the initial equations because all interatomic distances are at their equilibrium 
separation. From this equation, the momentum of atom number two can be calculated, 

thus allowing all other momenta to be calculated for their initial values via: 

(Pn+l ~ Pn) ^ 

=Eke 

2p 

In order to initiate a shock wave into the forty diatom chain, a ballistic particle 
(impact atom) with equal mass of a flourine atom imparts a potential upon the system. 
The interaction potential of the impact atom and with atom number one is given by the 
following equation: 



V- = 

* int 
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A = 0.0029 Hartree 
a=1.0(Bohrl) 

qi = position of chain atom number one 
Qo = position of impact atom 

The quantity A was chosen as twice the ground state energy value of a diatom pair, and 
a was chosen as 1.0 Bohr to optimize the trajectory of the impact atom in order to 
produce a movement of the atoms at a speed greater than 2500 m/sec (speed of a shock 
wave). The entire classical potential of the chain may now be written as a combination 
of the interaction potential from the ballistic particle, the sum of harmonic potential 
interaction of the particles and an anharmonic term; 



where 

qn = position of atom n 

Po = intramolecular equilibrium distance (Table 1) 

Ro = intermolecular equilibrium distance (Table 1) 

As the ballistic particle moves toward atom one with twice the initial momentum 
of atom number two, the repulsive potential becomes more significant, thus creating 
transverse motion within the chain. These atoms move according to Hamilton's 
equations of motion: 





40 



n=2,4,6 ... 
40 
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Upon numerical integration, these coupled differential equations yield atom momenta 
and positions as a function of time. The Adams-Moulton integrator was used in this 
study with the following parameters: 

Order = 6 

DT = 1 atomic time unit (2.5xl0'^ ps) 

EPS = 10'3, maximum relative degree of iteration 
N = 60,000 = numbers of time step iterations 
Results produce classical trajectory values for the one dimensional system. As stated 
previously, the two diatoms located at the center of the cluster are of main interest for 
this classical problem and for comparison to the quantum mechanical problem. 

In the quantum mechanical system, the four atoms are initially all located at their 
average equilibrium value corresponding to their counterparts in the classical system. 
These four atoms are connected to driving atoms (q< entering atom; q> absorbing atom) 
as formulated in the classical problem, and the motion of these driving atoms in the 
classical system is convened into the driving potential for the quantum system. This 
potential creates a perturbation within the quanutm model causing atom displacement. 
Since the one dimensional quantum mechanical system may be thought of as four 
harmonic or anharmonic coupled oscillators, motion in the system is described by four 
normal modes. The positions in the system may be expressed as a sum of initial 
cartesian position plus a time dependent displacement 5: 



qi =q° + 5i 

qi = Q2+ §2 

Q3 = q3 + 83 



(2c) 



(2b) 



(2a) 



™li 
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q4 = 94 + §4 (2d) 

From these relationships, calculations of transverse atom movement are performed. 

Formulation of the Hamiltonian for the Cluster 

For a set of atoms which are interconnected via sets of anharmonic springs, 
there are many types of internal motions, vibrations, which are dependent upon the 
initial displacements of the atoms. The Hamiltonian for the anharmonic system may be 
expressed as the sum of a harmonic part (H°), anharmonic part (Va), and time 
dependent part (V(t)). 

Ht = HO + Va + V(t) 

Formulation of Harmonic Hamiltonian 

The harmonic Hamiltonian for this one dimensional two diatom case is 
expressed as the sum of the potential and kinedc operators'^ 

H° = ^ [(q 1 - q< (0 - Ro)^ + (q> (0 - 94 ' Ro)^ 1 

2 2 
+ “ [(92* 9l ■ Pj + (94 ■ 93 ■ Po)T + "^(93 • 92" Rq) 

where 

Pi = mj Qi = momentum of atom i 
mi = mass of atom i 
qi = position of atom i 
Po = equilibrium intramolecular distance 



Ro = equilibrium intermolecular distance 
K< = Force constant of intermolecular bond 
K> = Force constant of intramolecular bond 

As stated previously, the positions of the four atoms as a function of time may 
be described as the sum of an initial position and a displacement position (5). The 
positions of the impact and absorbing atoms (q< and q>) may also be expressed as an 
initial position plus a displacement: 



The harmonic Hamiltonian is time independent, thus 5<(t) and 5>(t) are equal to zero 
until a permrbation is introduced. Utilizing the previous relations of initial positions of 
atoms, it is clear that: 



q<(t) = q< + 5<(t) 
q>(t) = q> + 6>(t) 



(3a) 



(3b) 



o o 



q2-qi = Po 



o o 



q4-q3-Po 

0 0^ 

q3-q2 = Ro 

o o 

q>-q4 = Ro 

o ® o 

qi-q<=Ro 



This enables simplification of the harmonic Hamiltonian: 



( 4 ) 
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* 

H = ( 5 1+54) +-^[(82-5 j) +(84-5.3) 



K. 






]+ (83-82)^ 



Formulation of Harmonic Time Independent Hamiltonian Matrix 

If an atom has an original position q° and throughout time moves to a different 
position q(t) the resulting vector describes the displacement (5) from its equilibrium 
value . This may be represented in mass scaled displacement coordinates (qj). 

qi = (mi)l/25j 

In matrix notation, 11)1 and S are defined as a row vector, and is a square matrix 
with the root values of atom masses on the diagonal. 

q = 



The kinetic energy may now be represented as 




1 I 



or in matrix notation 



- 1 ’ ‘ , 



where q '■ is the transpose of the row vector q. 



The potential energy of the system may also be represented in the same manner. 



We expand V in a Taylor series in the mass scaled displacements: 



j=i ' 3 ti ^ i j 



d V 



where 



V = potential energy 



Vo = V(5=0) = 0 




( Expanded about 5=0, the minimun of V ) 



Thus, the harmonic potential energy may be expressed as 
i j arii^j 

where Tj is the mass scaled coordinate row vector, is its transpose (a column 
vector), and K = force constant matrix which may be represented as follows: 



K -K> 0 0 

-K> K -K< 0 
0 -K< K -K> 
0 0 -K> K 



K< is the weak intramolecular force constant 



K> is the strong intermolecular force constant 
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K = K> + K< 

In order to analyze the system, an orthogonal transformation is required such 
that = A'K Define ^ = H]jA and = A^ir|^ where ^ is the normal coordinate row 

vector. The kinetic and potential matrices may now be rewritten with this substitution. 

t = ^C(a‘))((A)C0 = ^^ 

^=SM1/2A 6 = ^AtM-l/2 

V = (i)^(AkA)^ ' = (i)l|A(A‘KA)A'll ' 



The A^KA matrix may be redefined as the diagonal normal mode frequency square 
matrix (Q'^). The potential energy is then 

V = J?02^t 

From these relations the normal mode frequencies have been established (Table 
2). The displacements of the atoms from their respective equilibrium positions may be 
represented as linear combinations of the normal modes 



S = ^AtM-l/2 

where a representation of the A matrix follows: 



A = 



-0.5 -0.5 0.5 -0.5 
-0.5 -0.5 -0.5 0.5 
-0.5 0.5 -0.5 -0.5 
-0.5 0.5 0.5 0.5 



During the analysis of the two diatom case the values of the A matrix were all 
approximately + 0.5 or - 0 . 5 . This enabled the ^ operators to be presented as simplified 
combinations of the displacements of atoms from their respective equilibrium positions. 
(See Fig. 2 for description of modes.) 



Cl = - 0 . 5 Vm( 6 i + 82 + 83 + 84) 


( 5 a) 


C2 = + 0 . 5 Vm(- 8 i - 82 + 83 + 84) 


( 5 b) 


C3 = - 1 - 0 . 5 Vm(8i - 82 - 83 + 84) 


( 5 c) 


C4 = + 0 . 5 Vm(- 8 i + 82 - 83 + 84) 


( 5 d) 



Formulation of Anharmonic Time Independent Hamiltonian 

If one considers an arbitrary potential V(q) between atoms n and n +1 which has 
a minimum at q=qo, die potential function V(q) may be expanded in a Taylor series 
about qoi^d 

V(q) = a +( 3 (qn+l - qn - Qo)^ + Y(qn+1 - qn - qo)^ 

Gamma (y) represents the physical size of the anharmonicity (deviarion from harmonic 
potential well ) and is negative. If one were to make the assumption that the (qn+i - qn 
- qo)^ term is sufficiendy small compared to the first terms of the expansion, a 
harmonic oscillator would define the system. The (qn+i - qn - qo)^ interaction of the 
potential function cannot be ignored in real systems; thus, this term has been added to 
the total Hamiltonian for the model. In order to generate actual physical models, these 
cubic anharmonicities between the two diatoms and interconnecting sets of atoms were 



introduced. The anharmonicities for the two molecules in the cluster are 7(64 - 63)3 and 
'/{&l - 61)^. We are assuming that anharmoniciries exist only between diatom pairs, and 
that only the cubic anharmonicity is significant. 



From equations (5), the following normal mode operator equations were 
derived: 



Anharmonicities within the intramolecular bonds couple only modes three and 
four. If the intermolecular bond is anharmonic, all four modes would couple. 
Expansion of these terms resulted in the following anharmonic potential operator: 



Formulation of Time Dependent Hamiltonian (V(t)) 

In the time independent Hamiltonian, the impact and absorbing atoms were 
stationary; thus, the 5<(t) and 5>(t) were equal to zero. Now movement of these atoms 
must take place in order to propagate a shock through the cluster. The terms two of the 
harmonic Hamiltonian containing these groups were as follows: 



7(64 - 63)^ = (C4 “ C3 ) 

7(52-6i)3 = 7m‘'''(C4+C3) 



3 



3 




( 6 ) 





Expansion of these term yields: 



pi 




ir^ 




• 


II 






1 






w...- 




(7) 



2 



[5i2 - 25 i 5<(t) + 5<2(t) + 5>(t)2 - 25>(t)54+ 842] 



The rime independent terms, (5i2, 842 ), are formulated in the rime independent 
Harmonic Hamiltonian Eq. ( 4 ), thus equation ( 7 ) may be seperated into time 
dependent and rime independent terms. The time dependent terms will produce our 
driving potential for the atoms in the cluster. Expansion of 5i and 84 in terms of 
normal mode operators (eqns 5) ^'s yields: 

^ [8<(t) m"*^ (Cl + C2 - C3 + C4) +5>(t) m'^" (Ci - C2 - ^3 - ^4) + 5>(t)2 + 8<2(t)] 

Movement of the impact and absorbing atoms therefore generates a change in position 
which results in a potential being developed. 

Definition of Raising and Lowering Operators 

The operators which correspond to momentum and coordinates must satisfy the 
commutation relation. In matrix mechanics, the one dimensional harmonic oscillator 



Hamiltonian is as follows; 
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In order to proceed with normal mode analysis, the following raising (a^) and lowering 
operators (a) are introduced. 
a = Ctp^+PC 

* 

at= a*p^+ P C 

a*P is chosen as a purely imaginary value such that aP* = -a*P 



i 




P = 



g 



Substitution of these values yields the following: 



r Ir 2 ^2 2 hQ 

a a — ^p^ 4- g ^ ] + -;y- 



t „ hg 
aa=H- — 






(8a) 






So the Hamiltonian of the system may be described as 



u t tig 
H = a a + 



(8b) 



The Hamiltonian may now be applied to the Schrbdinger equation in order to evaluate 



the eigenstates. 

hIe> = e Ie> 
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(a a+^^)|E> = E|E> 

a ^a I E> = (E - ) I E> 

such that 

<E I a a I E> = (E - ^ )<E I E> 

<E I E> = 1 

<E|a a|E> = (E-^^) 

This will provide a harmonic basis for the problem. The states may be labeled as 
numbers with corresponding energies evaluated from the numbered states. 

The raising and lowering operators may now be defined through their action 
upon the eigenstates. 

a I n> = V nRQ |n-l> 



aM n> = V (n+i>hQ |n+l> 

In the model problem, the normal mode operators were defined from equations ( 8 ): 



C = CZ-fia . ) 



(10a) 




(10a) 



Using the above relations for a and a^ the cubic anharmonic operators may now be 
evaluated 



= (• 



2m; Q 



-)’^^l(n,)^^^ Ui- 1> + (n;+ l)'^^|ni+ 1>] 
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Ciln;> = ( [(n;) (nj- 1) I nj - 2> + (n;) I ri; > + (rij + 1) | rij > 

2rrii£2j 

+ (nj+ l)^^^ni + 2)^^^ lni + 2>] 



r"' I ^l/2. ,xl/2. o^l/2 I - .. ,1/2. ,, , ,1/2, 

Cilni> = ( ) [(nj) (tij-l) (nj-2) lni-3>+((n;) (rii-l)+(nj) (n,+ l) 



2rrijQj 



+ (n;)^^V(ni+ 1)^^^ !rii- 1> + ((ti; Kri; + (n^ +2)) 1 ti; + 1> 

+ (nj + 1 ) + 2 )^^^ (iii + 3) ! rij + 3>] 



The original basis set chosen for the problem consisted of the product of four 
normal modes ( I ni> 1 n 2 > 1 n 3 > 1 n 4 >). The normal mode operator acts only upon its 
corresponding basis state (i.e. ^ 4 ^ operated only upon 1 n 4 >). The normal mode 
operators were used to formulate the anharmonic and time dependent Hamiltonian 
matrix elements. In the total Hamiltonian matrix, the harmonic terms correspond to the 
diagonal elements and the anharmonic terms correspond to the symmetric off-diagonal 
elements. 



. I r I / n , 1/25 , 1 

<ni ! Ciln,> = (- — ) [(nj) ^ + (n. + 1) J 



1/2, 



1/2, 



2m i Qj 



The anharmonic Hamiltonian matrix was generated and diagonalized to obtain the 
appropriate eigenvalues and eigenstates ( Table 3 ). 

Z^HZ = S 

where Z = eigenvector matrix 
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= transpose of eigenvector matrix 
H = diagonal eigenvalue matrix 



Propagation of Time Dependent Hamiltonian 

The following passages formulate the time dependent perturbation introduced 
into the one-dimensional system.^l The basic solution to the time dependent 
Schrodinger equation is as follows: 

Hvj/= (ih)(^) 

where H is a function of position and time in a Cartesian coordinate system. This is 
extremely difficult to solve; therefore, the Hamiltonian is divided into a time- 
independent and time dependent part: 

H = Hti + V(t) 

Hti = time independent tenn of Hamiltonian 
V(t) = time dependent term of Hamiltonian 
'P = \|/(t)e-i(E/fi)i 

The spatial equation to solve is 
Hxi\j/ = E\)/ 

E = separation constant 

It has well behaved solutions when E is equal to En 
HnVn® = EnVl/n° 

Assume the 'P(q,t) can be written as a linear combination of orthogonal functions: 

4^(q,t) = Xan(t)vI(q.t) 

n 

vj/ = vp 
^ n ^ n ^ 



The quantity is described by the eigenfunctions of a harmonic oscillator and the 

quantities an(t) are defined as the occupation coefficients. Since the wave function is 
normalized, 

£ a„(t)X(t)=l 

n 



the series expansion for T(q,t) may now be substituted into the Schrodinger wave 
equation yielding; 







where m = 1, 2, 3, ... 

This yields the rate of change (with respect to time) of the m^ component of the 
basis. Am (rate of amplitude change) depends upon all basis amplitudes. The 
following equations are shown to exemplify the above differential equation. 



h dai 
i dt 



= 3i 



/ 



-h da; 
i dt 




V(t) dt + a2 



O* O 

V(t) dt + a2 



/ 




O* 



o 

V(t)\|f2dt + ... + aK 



O 

V(t)\|/2dt + ... + aK 



/ 

/ 



O * O 

Vi V(t)\[/^dt + ... 



o* o 

V 2 V(t)\|/f.dt+... 



UtiUzing these coupled first order differential equations, numerical analysis using the 
Adams Moulton integrator provided the coefficients for the various states of the model. 



Calculation of Quantum Position of Atoms 

In order to calculate positions of atoms, computations of expectation values of 
displacements from equilibrium were performed. The matrix elements of operators 5j 
were evaluated in the normal mode basis, and then transformed into the harmonic 
representation using the eigenvector matrices , as follows: 

Z^DjZ = Dj 

where 

Z^ = transpose of eigenvector matrix 

Z = eigenvector matrix 

Dj = displacement matrix 

D'i = transformed displacement matrix 

Upon transformation, D'i was then summed over each state with its coiresponding 
energy and probability: 

?‘(t)D'i ?(t) 

where 

[P(t)]i = ci(t) exp [-iEit/Ii] 

where P is a column vector based on probablity and eigenvalues, and is a row 
vector. This summation yields the displacement of atom i from its equilibrium 
position. 

Calculation of Quantum Bond Energies 

The total energy of a particular bond is expressed as the summation of the 
bond's potential and kinetic energy. The potential energy for the quantum system is 
calculated from the expectation values of position for the corresponding atoms. 
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Potential energy = (l/ 2 )(Kbond)(Sj - 5i)^ 

From equations (5) the potential energy can be represented in terms of normal modes: 

PE, 2 = ij) (K^ ( 52 - 5 ,) 2 = ( 1 ) (K^ 

PE23 = (^) (K^ (83 - 82)2 = ( 1 ) (K<) (C2 - 
PE34 = (^) (K>) (84 - 83)2 = ( 1 ) (K^ (C4 + C3)" 

These operators will act upon our normal mode basis set in order to establish potential 
energy matrices. Upon transformation of the potential energ>' matrix with the 
eigenvector matrix, summation over all states at a particular time yields the potential 
energy of a bond. 

Z^EjjZ = PE,j 

? \t) PEjj P(t) = potential energy of bond ij 
As stated previously the kinetic energy is calculated by: 

- 1 ‘ ‘ , 

where C * is a column vector and is normal mode momentum. As mentioned with 
respect to position, normal mode momentum ( P^i ) may be expressed as a summation 
of corresponding atom momentum (pi) 

p^j = -0.5Vm (Pi + P2 + P3 + P 4 ) = (l/(i^)) (ai-a 1 ) 

P(^2= (-pi - P2 + P3 +P 4 ) = ( l/(iV^)) (^ 2 -^ 2 ) 
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P ;3 = + 0 . 5 Vm (Pi - P 2 - P 3 + P4) = (l/(il^)) (33-33) 
p^4 = -t- 0 . 5 Vm (-pi + P2- P 3 + P4) = (l/(if 2 )) (a4-a4) 

The kinetic energy of the i-j bond is expressed as: 



T 

‘j o 

2|l 



which is expressed by normal momentum operators as: 



(PC4 ' PC3 ) 

T 1 2 = = kinetic energy operator for bond R 1 2 

2fi 



(P;2‘ P^4)' 



T23 = = kinetic energy operator for bond R23 



2(1 



(PC 4 PC3 ) 

T34 = — ^ ^ = kineric energy operator for bond R34 

2(t 



Each kinetic energy operator established a kinetic energy matrix which is transformed to 
the harmonic representation via the eigenvector matrix Z. 

z‘KS;jZ = KE;j 

Upon transformation we sum over all states: 

p\t)KEijP(t) 
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Calculation of Classical Bond Energies 

Classical potential and kinetic bond energies are calculated directly from atom 
position. 

Potential energy = (f)(Kbond)(qn+l - qn - po)^ 

Kinetic energy = C/2^) (Pn+i ' P) 

Calculation of Atomic Period 

The atomic period for this model was based upon the stonger force constant. 
C0 = (K>/p)1/2 
T = 27t / 0 ) = 2. 1 X 10 3 atu 



Energy Evaluations of the Two Diatom Model 

In order to evaluate the normal mode energies, the original Hamiltonian matrix 
was separated into four normal modes, a time-independent anharmonic term, and the 

time-dependent potential term. 

nj=I9Ij-i-Pii2 + Ki2 + ]KI^ + Vj|+V(t) 

Each matrix was initially transformed into the harmonic basis by multiplication with the 
eigenvector and transpose eigenvector matrix (Z = eigenvector matrix). 

Z^HjZ =H; 

Upon transformation, each term was summed over all possible states to provide the 
energies of the particular parts. 

P \t) P3 i P(t) 



III. Discussion of Results 



A. Classical Trajectories 

When the ballistic particle imposes a shock, it adds .005884 Hartrees of energy 
to the chain of atoms. The shock wave propagates through the classical chain at 
approximately 6.5 x lO'^ ao/ps, causing displacement of atoms ( Figures 3,4). The 
classical system responds as expected with respect to atom position and bond energies. 
Atom position is directly related to the neighbor atoms, and atoms tend to oscillate in 
diatom pairs. The impact atom provides a disturbance within the lattice which is not felt 
by the central atom pairs for approximately eight periods (Figs. 5). The transverse 
displacement of the entering atom (qo) imparts energy to the first atom in the lattice 
which causes the transverse movement of this atom and so fonh; thus, resulting in the 
propagation of the disturbance throughout the lattice. Since this shock wave is a 
comparatively low energy disturbance with respect to previous calculations 
[Wyatt,Marston], atoms are not significantly displaced and atom separation is fairly 
constant. The atoms connected by a strong bond tend to move as a pair, with most of 
the lattice deformation occuring within the intermolecular bond. 

Bond energies indicate that energy is transferred between kinetic and potential 
energy . As the energy from the shock wave is transferred to the bonds, 
nonuniformities in the oscillator behavior are observed. The sffong bonds exhibit a 
tendency to remain relatively constant, while most of the energy from the shock is 
transferred to the intermolecular bond ( Figures 6-11 ). 
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B . Convergence of the Quantum Model 

The time dependent coefficients for the quantum states aj (t) will help us in 
examining the convergence of the results of our system, the nature of occupation of 
various quantum states, and the conservation of probability. Quantum eigenvalues of a 
sixteen state model are listed in Table 3. In the present model, the classical positions of 
q< and q> are used in calculating a driving potential. The classical system was 
initialized with all atoms in the chain located at their respective eqilibrium values. Initial 
momenta were calculated for these atoms, and Hamilton's equations of motion were 
solved for position and momentum. Upon solving the equations of motion, the atoms 
were no longer located at equilibrium displacements. This has been defined as t = 0 for 
the system. Since q< and q> were displaced from equilibria, a small driving potential 
for the system was established causing population of excited states in the quantum 
model. The model may now be thought of as being non-zero temperature. ( atoms are 
oscillatory ). 

When a shock wave is introduced to the quantum model, energy is transferred 
into the crystal lattice resulting in excitation and population of higher energy states. 

This corresponding population of states may be represented as ai*ai where ai is the 
coefficient of occupation of state i ( .04% deviation in the 16 state model, .008% 
deviation in the 100 state model). Plots of a 16 state basis illustrate a periodic transfer 
of energy from state to state and a tendency to retain most of the energy within the 
lower mode states (Figures 12-27 ). Harmonic and anharmonic 16 state bases exhibit 
similar tendencies. Examination of a larger basis (100 states) shows a larger retention 
of energy in lower mode states. Since there is a ratio of 5:5:2:2 (ni:n 2 :n 3 :ri 4 ) in the 100 
state system compared to a 2:2:2:2 ratio in our lower (16 state) model, the energy 
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distributes itself differently; thus, allowing for even less population of higher energy 
levels (Figures 28-33). 

C . Comparison of Position Trajectories 

Comparison of expectation values for the positions of atoms for several bases 
indicates fairly good convergence of the system. (Figures 34-41 ) We have learned 
from [Wyatt and Marston] that quantum harmonic atom positions at high energy 
converges towards the classical results. Although there are deviations from this 
expected result one can see that the larger harmonic basis sets move toward 
convergence to classical results. (Figures 42-45 ) This is also the case with 
convergence of the anharmonic system. 

Anharmonicities produce deviations in energy level spacing as compared with 
a purely harmonic system. The deviations in energy spacing produces a variance in 
energy state occupation; thus, there will be variances in position plots of harmonic and 
anharmonic models. Since there are only two states in the higher frequency modes, the 
100 state system does not fully describe the cubic anharmonicity making differences 
between the harmonic and anharmonic models small. (Figures 46-49 ) Atom motion in 
the anharmonic model deviates slightly in frequency and amplitude as compared to the 
harmonic model. 

In comparison to classical harmonic trajectories, the quantum harmonic model 
does not converge exactly to the classical harmonic trajectory, as expected [Wyatt, 
Marston]. The initial oscillatory motion of the quantum model prior to a shock is 
deviates slightly coupled with small deviations of average position away from average 
equilibrium values. This result is most likely caused by the low number of energy 
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states in the quantum mtxiel. If one were to examine the frequencies of the four modes 
(Table 2), one can see that they are in a ratio of 11.5:6.6:1:1 

The quantum model currendy has a 5:5:2:2 basis ratio. Movement in each 
system should be characterized by oscillations of diatom pairs. For instance, normal 
oscillations of atoms one and two should be directly out of phase as observed in the 
classical system; however, careful observation of quantum diatom pairs yields small 
phase dissimilarities which may be caused by the population of various energy states. 
Since we have a low number of internal mode states (modes one and two), the higher 
modes are increasingly populated. If one were to examine the modes (Figure 2), one 
can see that linear combinations of the modes describe the normal oscillator motion of 
the atoms. Since the higher energy modes are more populated, we observe phase 
dissimilarities within diatom pairs. Classical anharmonic results versus quantum 
anharmonic results correspond in a similar fashion. 

The quantum harmonic model of [Wyatt, Marston] was easily separated into 
four independent modes. This allowed for vast numbers of low frequency mode states 
as compared to the higher energy modes (500:200:100:100) (ni:n 2 :n 3 :ri 4 ) giving a total 
basis size of 10^. The low frequncy modes were now capable of describing center of 
mass and oscillatory motion. Addition of this number of states to the present model 
makes the problem numerically intractible. 

It is also interesting to note the dependence of the present model upon q>. The 
absence of the potential created by q> causes great phase dissimilarities with respect to 
diatom oscillations, and it also inhibits proper atom atom displacements with respect to 
a shock induced lattice. Several cases have been examined ( Figures 50-59 ) in order to 
show these effects. Since the q> potential is no longer present prior to a shock ( ramp 
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shock ) or not present at all, this causes the excited state population to vary 
significandy; thus, atom energy and displacement become less uniform. The potential 
caused by q> has a significant affect upon lattice dynamics. 

Another interesting point to note is the dependence of this model on the size of 
the anharmonicity. The degree of anharmonicity was increased 10 fold as a trial 
experiment Small deviations in position are observed as well as small frequency shifts 
in oscillations Fig. ( 60-61 ). In the present model, we are limited to a ground state and 
first excited state in modes three and four. This limitation in basis size effectively 
reduced the survival terms of the anharmonic operators; thus, allowing for small 
deviations in frequency and amplitude of atom positions. 

D. Quantum Mode Energy Comparisons 

A slightly larger difference between our quanutm harmonic and anharmonic 
models is observed in mode energies. (Figures 62-65 ). Energy level spacing in the 
harmonic model is different as compared with the anharmonic model in which energy 
level spacing decreases towards higher levels. This spacing difference exhibits 
variances in mode energies prior to and after shock wave introduction. If both models 
were initially stationary states, one would observe that mode energies would start at 0.5 
quanta. As stated earlier, t =0 corresponds to a time in which atom positions are 
displaced from their equilibrium values causing a driving potential to be created by q< 
and q>. This deviation in q< and q> is fed into our quantum system providing for 
population of some excited states. Since the excited states became populated, initial 
values of mode energies increased above 0.5 quanta. Initially, the harmonic model 
distributes energy to its excited states in a different fashion as compared with the 



anharmonic model. The anharmonic model slowly distributes its energy to higher 
states then reaches a steady value prior to the shock wave. Energy transitions are 
much smoother in the anharmonic model because energy levels have been lowered. 

As shock energy is introduced to both systems, the anharmonic model reacts 
more quickly and distributes energy to higher states. This is evidenced by the fact that 
the higher modes reach a steady state at slightly higher values as compared to the 
harmonic model, and the lower mode energies in the anharmonic case achieve smaller 
energies as compared to the lower modes. In each case, the upper modes pick up 
approximately 0.2 quanta of energy as compared with approximately 0.6 quanta for the 
lower modes. A ratio of this (0.6;0.2) corresponds quite well with a ratio of states 
(5:2). 

Addition of states to the lower frequency modes in these models will definitely 
improve the energy transfer prior to and during the shock wave. Energy will be 
distribute differently among the lower states allowing for the system to stabilize. 

E. Comparison of Bond Energies 

The classical models react as expected with respect to bond energies. Diatom 
pairs are oscillating in a uniform fashion prior to a shock wave resulting in kinetic and 
potential energy to be directly out of phase and total energy to remain constant. As the 
system is perturbed, oscillations and dispacements vary; thus, variances in potential, 
kinetic and total energy may be observed. The strong bonds have an inherent tendency 
to resist change with respect to a low energy shock. ( Figures 6-11 ) 
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Prior to a shock in the quantum models, energy within the bonds is dependent 
upon the coupling of various modes. Since the system is non-stationary, population of 
various energy levels vary with time. This causes individual bond energies to fluctuate 
transferring energy between bonds. As the shock wave is introduced to the system, 
bonds exchange energy slightly faster and the overall energy within the bonds is higher 
( Figures 70-79 ). There are few differences between the harmonic and anharmonic 
quantum models; however total bond energy is slighdy lower in the anharmonic 
system. Careful examination of the potential versus kinetic energy plots of bond 
energy in both quantum systems shows that potential and kinetic energies are in phase. 
Part of this is due to the sharing of energy between bonds; however, some is caused by 
an small number of low frequency mode states (modes one and two ). Excessive 
energy is placed into the higher mode states; thus causing an imbalance in the system. 
This was not as evident in positions as it is in bond energies. 



IV. Conclusions 



Examination of a crystal lattice using a coupled mode basis yielded a great deal 
of information concerning both harmonic and anharmonic models. The following 
conclusions may be drawn from these results; 

(1) The deviations in energy spacing from the anhaimonicities cause 
transitions to occur more easily as compared with the harmonic model 

(2) Addition of states to the quantum model allowed for convergence toward 
classical results. 

(3) The anharmonicities produce slight deviations in frequency of oscillation 
and displacement because of energy spacing differences. 

(4) The quantum model is not only dependent on the potential created by q< ; 
the potential created by q> is necessary to allow the atoms to act as a 
cluster. Any deviations in the potential created by q> causes the model to 
change significantly. 

(5) Energy spacing differences allowed for deviations in mode energies. 

(6) Bond energies within the classical models correspond directly to atom 
positions. There are instances in time in which potential energy is at a 
maximum (corresponding to maximum deviation of atom position from 
equilibria) and minimum kinetic energy. Total energy is constant prior to 
a shock. Most of the energy is transferred to the intermolecular bond. 

(7) Bond energies within both quantum models are not constant prior to a 
shock. Since energies are dependent upon mode population, total energy 
within a bond may oscillate. Energy may be transfered from bond to 
bond. The total energy of all bonds prior to a shock should be constant. 
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(8) Deviations from the expected results of [Wyatt, Marston] in a harmonic 
model were caused primarily by the absence of states which resulted in 
normal oscillator motion of the atoms and center of mass motion to be 
inconsistent. Careful examination of the mode frequencies yielded 
information concerning the proper number of states one must have in 
order to describe atom motions properly. Lack of a substantial number of 
states also restricted the amount of energy transferred into the lattice. 

Upon re-examination of the problem, one sees that the total quantum 
Hamiltonian was written as follows; 

Hj = HO +Va +V(t) 

The anharmonic substituent (Va) of the Hamiltonian was formulated using a coupled 
normal mode operator: 

V, = 7(m)-^'^2C4+6C3C4) 

Instead of using a Hamiltonian which had all four modes coupled, we should have 
formulated the problem using a harmonic uncoupled Hamiltonian for modes one and 
two, and a coupled anharmonic Hamiltonian for modes three and four. 

Rj= Hi + H2+ H34+ V34+ Vi(t) + V2(t) + V34(t) 

where 

Hj= total Hamiltonian 

H° = Harmonic Hamiltonian for mode i 

H34= Harmonic Hamiltonian of coupled modes three and four 

V34= Anharmonic Hamiltonian of coupled modes three and four 

Vj (t) = Potential for mode i 

V34 (t) = Potential for coupled modes three and four 



This would have allowed for an increased number of low frequency mode state (modes 
one and two); thus, placing the energy levels of all the modes in proper ratios. 

In conclusion, an entirely coupled Hamiltonian severely restricted proper 
formulation of this problem; however, reformulation using uncoupled modes should 
prove beneficial. Since desired results were not obtained using the harmonic model, 
deviations are most likely present in our anharmonic model. Once the Hamiltonian is 
uncoupled, an extremely large number of states may be added, and it will afford 
extremely large additions of energy to the crystal lattice. 
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Table 1 

Equilibrium values for position, atom separations, mass, 
and force constants 



q<= -6.9316 a„ 
q° = -3.4658 a^, 
q^= -1.7329 a, 

p„ = +1.7329 a, 

mi = m 3 = 34629 amu 
m 2 = m 4 = 34629 amu 



q° = +6.9316 a^ 
q 3 = +1.7329 a„ 
q° = +3.4658 a„ 

R„ = +3.4658 a„ 

K< = 0.00468 Hartree/Bohr^ 
K> = 0.14985 Hartree/Bohr^ 



|l = 17314.5 amu 



= -3.0 X 10"^ Hartree/Bohr^ 



Table 2 



Frequencies and energy spacings for the normal modes. 



Normal Mode Frequencies (a.u.) 


Energy Spacing (hf2)(eV) 


Ql = 0.000256 


6.966 X 10-3 


^2 = 0.000445 


12.108 X 10-3 


^3 = 0.002927 


79.644 X 10-3 



Ot = 0.002949 



80.243 X 10-3 
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Table 3 

Eigenvalues for the anharmonic and harmonic quantum 16 state system. 
State (na 03 02 ni) Anharmonic Energy (a.u.) Harmonic Energy (a.u.) 



0000 


0.003284 


0.003290 


0001 


0.003540 


0.003546 


0010 


0.003722 


0.003735 


0011 


0.003989 


0.003992 


0100 


0.006163 


0.006217 


1000 


0.006235 


0.006239 


0101 


0.006435 


0.006473 


1001 


0.006594 


0.006496 


0110 


0.006634 


0.006662 


1010 


0.006684 


0.006685 


0111 


0.006896 


0.006919 


1011 


0.006942 


0.006942 


1100 


0.009165 


0.009166 


1101 


0.009421 


0.009423 


1110 


0.009610 


0.009612 


nil 


0.009859 


0.009869 
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Figure 1 



Figure 2 

Figure 3 

Figure 4 

Figure 5 
Figure 6 
Figure 7 
Figure 8 



Figure Captions 



Representation of the four atom cluster , with q< and q>. 

q< = initial cartisian coordinate position of entering qroup atom 

qi = initial cartesian coordinate position of atom one 

<\2 - initial cartesian coordinate position of atom two 

q 3 = initial cartesian coordinate position of atom three 

qa = initial cartesian coordinate position of atom four 

q> = initial cartesian coordinate position of the absorbing atom 

K> = force constant utilized between diatom pairs 

K< = force constant utilized for entering, absorbing, and coupling of 

diatoms 

Rp = cartesian distance utilized for entering, leaving, and coupling of 
diatoms 

Po = cartesian distance utilized for coupling of diatom pairs 

Pictoral representation of the four normal modes. 

Normal mode one represents center of mass motion of the four atoms. 
Normal mode two represents paired diatom stretching about the central 
bond between atoms two and three. 

Normal mode three repesents an attractive stretching between atoms one 
and two, and a repulsive interaction between atoms three and four. 
Normal mode four represents a repulsive interaction between atoms one 
and two, and an attractive interaction between atoms two and three. 



Classical harmonic position of atoms q<, qi, q2, q3, qa, and q> . 
Positions are measured in atomic units, periods correspond to one 
harmonic oscillation of a diatom pair. Atom positions are measured 
from the center of mass for the chain. 

Classical anharmonic position of atoms q<, qi, q2, q3, qa, and q> . 
Positions are measured in atomic units, periods correspond to one 
harmonic oscillation of a diatom pair. Atom positions are measured 
from the center of mass for the chain. 

Position of the ballistic particle. Atom position is measured in atomic 
units. Position is measured from the center of mass for the chain. 

Classical harmonic bond energy of the intramolecular bond ( R12 ) 
between atoms qi and q2. Energy is measured in harmonic quanta. 

Classical harmonic bond energy of the intramolecular bond ( R34 ) 
between atoms q 3 and qa. Energy is measured in harmonic quanta. 

Classical harmonic bond energy of the intermolecular bond ( R23 ) 
between atoms q2 and q 3 . Energy is measured in harmonic quanta. 
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Figure 9 
Figure 10 
Figure 1 1 
Figure 12 

Figure 13 

Figure 14 

Figure 15 

Figure 16 

Figure 17 

Figure 18 

Figure 19 

Figure 20 

Figure 21 



Classical anharmonic bond energy of the intramolecular bond ( R12 ) 
between atoms qi and qi- Energy is measured in harmonic quanta. 

Classical anharmonic bond energy of the intramolecular bond ( R34 ) 
between atoms q3 and q4. Energy is measured in harmonic quanta. 

Classical anharmonic bond energy of the intermolecular bond ( R23 ) 
between atoms q2 and q3. Energy is measured in harmonic quanta. 

Representation of the occupation of the two lowest energy states in the 
16 state quantum harmonic model. The probablities (nq n3 n2 ni) 
correspond to the quantum number of a particular mode. 

Representation of the occupation of energy states (0010) and (0011) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (U4 n3 n2 n|) of a panicular mode. 

Representation of the occupation of energy states (0100) and (1000) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (r4 n3 n2 ni) of a particular mode. 

Representation of the occupation of energy states (0101) and (1001) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (114 n3 n2 ni) of a particular mode. 

Representation of the occupation of energy states (0110) and (1010) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (U4 U3 n2 ni) of a particular mode. 

Representation of the occupation of energy states (0111) and (101 1) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (nq U3 U2 ni) of a particular mode. 

Representation of the occupation of energy states (1 100) and (1 101) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (04 n3 U2 ni) of a particular mode. 

Representation of the occupation of energy states (1110) and (1 1 1 1) in 
the 16 state quantum harmonic model. The numbers correspond to 
quantum numbers (ri4 n3 n2 ni) of a particular mode. 

Representation of the occupation of the two lowest energy states (0000) 
and (0001) in the 16 state quantum anharmonic model. TTie numbers 
correspond to quantum numbers (n4 n3 n2 ni) of a particular mode. 

Representation of the occupation of energy states (0010) and (0011) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (04 n3 U2 ni) of a particular mode. 
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Figure 22 



Figure 23 



Figure 24 



Figure 25 



Figure 26 



Figure 27 



Figure 28 



Figure 29 



Figure 30 



Figure 31 



Figure 32 



Figure 33 



Figure 34 



Representation of the occupation of energy states (0100) and (1000) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (ri 4 U 3 U 2 nj) of a particular mode. 

Representation of the occupation of energy states (0101) and (1001) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (04 U 3 U 2 nj) of a particular mode. 

Representation of the occupation of energy states (0110) and (1010) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n 4 U 3 U 2 ni) of a particular mode. 

Representation of the occupation of energy states (0111) and (101 1) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (n 4 n 3 n 2 nj) of a particular mode. 

Representation of the occupation of energy states (1 100) and (1 101) in 
the 16 state quantum anharmonic m.odel. The numbers correspond to 
quantum numbers (04 U 3 U 2 nj) of a particular mode. 

Representation of the occupation of energy states (1 1 10) and (1111) in 
the 16 state quantum anharmonic model. The numbers correspond to 
quantum numbers (04 U 3 U 2 ni) of a particular mode. 

Representation of the ground state occupation of (0000) in the 100 state 
quantum anharmonic model. The numbers correspond to quantum 
numbers (rt 4 n 3 n 2 ni) of a particular mode. 

Representation of the occupation of the first e.xcited state (0001) in the 
100 state quantum anharmonic model. The numbers correspond to 
quantum numbers (a 4 n 3 n 2 ni) of a particular mode. 

Representation of the occupation of excited energy state (0010) in the 
100 state quantum anharmonic model. The numters correspond to 
quantum numbers (ri 4 n 3 U 2 ni) of a particular mode. 

Representation of the occupation of excited energy state (1134) in the 
100 state quantum anharmonic model. The numters correspond to 
quantum numbers (U 4 U 3 U 2 ni) of a particular mode. 

Representation of the occupation of excited energy state (1143) in the 
100 state quantum anharmonic model. The numters correspond to 
quantum numbers (04 U 3 U 2 nj) of a particular mode. 

Representation of the occupation of excited energy state (1 144) in the 
100 state quantum anharmonic model. The numters correspond to 
quantum numbers (U 4 U 3 n 2 ni) of a particular mode. 

Comparison of the positions of atoms one and two (qi,q2) in the 16 
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Figure 35 



Figure 36 



Figure 37 



Figure 38 



Figure 39 



Figure 40 



Figure 4 1 



Figure 42 



Figure 43 



state and 100 state quantum harmonic models. Postions are measured in 
atomic units (ao) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 16 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (Uq) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms one and two (qi,q2) in the 64 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (a<}) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 64 
state and 100 state quantum harmonic models. Postions are measured in 
atomic units (ao) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms one and two (qi,q2) in the 16 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (uq) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 16 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (Uq) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms one and two (qi,q2) in the 64 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (Uq) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 64 
state and 100 state quantum anharmonic models. Postions are measured 
in atomic units (ao) and one period corresponds to one oscillation of an 
atom. 

Comparison of the positions of atoms one and two (qi,q2) in the 100 
state quantum harmonic model and the classical harmonic model. 
Postions are measured in atomic units (ao) and one period corresponds 
to one oscillation of an atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 100 
state quantum harmonic model and the classical harmonic model. 
Postions are measured in atomic units (uq) and one period corresponds 
to one oscillation of an atom. 
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Figure 44 

Figure 45 

Figure 46 
Figure 47 
Figure 48 
Figure 49 

Figure 50 

Figure 51 

Figure 52 

Figure 53 

Figure 54 

Figure 55 

Figure 56 



Comparison of the positions of atoms one and two (qi,q2) in the 100 
state quantum anharmonic model and the classical anharmonic model. 
Postions are measured in atomic units (ao) and one period corresponds 
to one oscillation of an atom. 

Comparison of the positions of atoms three and four (q3,q4) in the 100 
state quantum anharmonic model and the classical anharmonic model. 
Postions are measured in atomic units (ao) and one period corresponds 
to one oscillation of an atom. 

Varations in atomic position (ao) of atom one (qi) in the quantum 100 
state anharmonic and harmonic models. 

Varations in atomic position (ao) of atom two (q2) in the quantum 100 
state anharmonic and harmonic models. 

Varations in atomic position (ao) of atom three (q3) in the quantum 100 
state anharmonic and harmonic models. 

Varations in atomic position (ao) of atom four (q4) in the quantum 100 
state anharmonic and harmonic models. 



Positions of atoms qi and q2 in the 16 state harmonic model in which 
the potential created by q> is zero at all times. Positions are measured in 
atomic units. 

Positions of atoms q3 and q4 in the 16 state harmonic model in which 
the potential created by q> is zero at all times. Positions are measured in 
atomic units. 

Positions of atoms qi and q2 in the 16 state harmonic model ; the 
potential created by q> is applied as a step jump at t = 10 periods. 
Positions are in atomic units. 

Positions of atoms q3 and q4 in the 16 state harmonic model ; the 
potential created by q> is applied as a step jump at t = 10 period. 
Positions are measured in atomic units. 

Positions of atoms qi and q2 in the 16 state harmonic model ; the 
potential created by q> is applied as a ramp increase between t = 10 and 
12.5 periods. Positions are measured in atomic units. 

Positions of atoms q3 and qa in the 16 state harmonic model ; the 
potential created by q> is applied as a ramp increase between t = 10 and 
12.5 periods. Positions are measured in atomic units. 

Positions of qi and q2 in the 16 state harmonic model ; the potential 
created by q> is the actual potential of q> shifted 3 periods earlier. 
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Figure 57 
Figure 58 
Figure 59 
Figure 60 

Figure 6 1 

Figure 62 
Figure 63 
Figure 64 
Figure 65 
Figure 66 
Figure 67 
Figure 68 
Figure 69 
Figure 70 



Positions are measured in atomic units. 

Positions of q 3 and q 4 in the 16 state harmpnic model ; the potential 
created by q> is the actual potential of q> shifted 3 periods earlier. 
Positions are measured in atomic units. 

Positions of qi and q 2 in the 16 state harmonic model ; the potential 
created by q> is the actual potential of q> shifted 3 periods later. 
Positions are measured in atomic units. 

Positions of q 3 and q 4 in the 16 state harmonic model ; the potential 
created by q> is the actual potential of q> shifted 3 periods later. 
Positions are measured in atomic units. 

Positions of atoms qi and q 2 in the 100 state quantum anharmonic 
models. The solid line plots are atom positions with gamma = -3.0 x 10 
-4 Hartree/ Bohr 3. The dashed line plots indicate atom positions with 
ten times the anharmonicity. 

Positions of atoms q 3 and q 4 in the 100 state quantum anharmonic 
models. The solid line plots are atom positions with gamma = -3.0 x 10 
-4 Hartree/ Bohr 3. The dashed line plots indicate atom positions with 
ten times the anharmonicity. 

Distribution of energy into mode one of the 100 state quantum harmonic 
model. Energy is measured in quanta. 

Distribution of energy into mode two of the 100 state quantum harmonic 
model. Energy is measured in quanta. 

Distribution of energy into mode three of the 100 state quantum 
harmonic model. Energy is measured in quanta. 

Distribution of energy into mode four of the 100 state quantum 
harmonic model. Energy is measured in quanta. 

Distribution of energy into mode one of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 

Distribution of energy into mode two of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 

Distribution of energy into mode three of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 

Distribution of energy into mode four of the 100 state quantum 
anharmonic model. Energy is measured in quanta. 

Potential and kinetic energy of bond ( R12 ) between atoms qi and q 2 in 
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Figure 7 1 
Figure 72 
Figure 73 
Figure 74 
Figure 75 

Figure 76 

Figure 77 
Figure 78 

Figure 79 



the 100 state quantum harmonic model. Energy is measured in quanta. 

Potential and kinetic energy of bond ( R34 ) between atoms q3 and q4 in 
the 100 state quantum harmonic model. Energy is measured in quanta. 

Total energy of bonds ( R12 ) and ( R34 ) in the 100 state quantum 
harmonic model. Energy is measured in quanta. 

Potential and kinetic energy of bond ( R23 ) between atoms q2 and q3 in 
the 100 state quantum harmonic model. Energy is measure in quanta. 

Total energy of bond ( R23 ) between atoms q2 and q3 in the 100 state 
quantum harmonic model. Energy is measured in quanta. 

Potential and kinetic energy in bond ( R12 ) between atoms qi and q2 in 
the 100 state quantum anharmonic model. Energy is measured in 
quanta. 

Potential and kinetic energy in bond ( R34 ) between atoms q3 and qa in 
the 100 state quantum anharmonic model. Energy is measured in 
quanta. 

Total energy in bonds ( R12 ) and ( R34 ) in the 100 state quantum 
anharmonic model. Energy is measured in quanta. 

Potential and kinetic energy of bond ( R23 ) between atoms q2 and q3 in 
the 100 state quantum anaharmonic model. Energy is measured in 
quanta. 

Total energy in bond ( R23 ) between atoms q2 and q3 in the 100 state 
quantum anharmonic model. Energy is measured in quanta. 
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